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ABSTRACT 

Recent observations with the Extreme Ultraviolet Explorer (EUVE) suggest that at least some clusters 
of galaxies are luminous sources of extreme ultraviolet (EUV) radiation. It is not clear yet whether EUV 
emission is a general feature of clusters; for the purposes of limiting the contribution to the background 
radiation, we assume that it is true of most clusters. Assuming that the source of the EUV emission is 
inverse Compton (IC) scattering of the Cosmic Microwave Background photons by relativistic electrons, 
we construct a simple model for the expected average emission from clusters as a function of their mass 
and the redshift of interest. Press-Schechter theory is used to determine the abundance of clusters of 
various masses as a function of redshift. We determine the amount of background radiation produced by 
clusters. The total mean intensity, spectrum, and the ionization rates for H I and He II are determined 
at present and at a variety of redshifts. Because clusters form by the merger of smaller subclusters, the 
amount of EUV background radiation should be larger at present than in the past. We compare our 
results to the ionizing background expected from quasars. We find that while clusters do contribute a 
significant EUV background, it is less than a percent of that expected from quasars. 

Subject headings: cosmic rays — diffuse radiation — galaxies: clusters: general — intergalactic 
medium — large-scale structure of the universe — ultraviolet: general 



1. INTRODUCTION 

Recently, extreme ultraviolet (EUV) and very soft X- 
ray emission has been detected from a number of clus- 
ters of galaxies. Except for the Coma and Virgo cluster 
(Bowyer, Lampton, & Lieu 1996; Lieu et al. 1996a,b,c; 
Bowycr, Berghofer, & Korpela 1999; Berghofer, Bowyer, 
& Korpela 2000), the individual detections are controver- 
sial. Other clusters which may have been detected in- 
clude Abell 1795 and Abell 2199 (Bowyer, Lieu, & Mittaz 
1998; Mittaz, Lieu, & Lockman 1998), Abell 4038 (Bowyer 
et al. 1998), and Abell 4059 (Bowyer 2000). The excess 
EUV detections in the rich clusters (i.e., not Virgo) corre- 
spond to luminosities of ~ 10 44 ergs s _1 . The detections 
of Abell 1795 and Abell 2199 as extended EUV sources 
have been quite controversial. For both clusters, detec- 
tions with EUVE were claimed by Bowyer et al. (1998) 
and Mittaz et al. (1998). However, Bowyer, Berghofer, 
& Korpela (1999) argued these detections were due to an 
incorrect subtraction of the background. In Abell 2199, 
Lieu et al. (1999b) reobserved the cluster, and continue 
to claim detection of an extended EUV source. An associ- 
ated soft X-ray excess may have been seen with BeppoSAX 
(Kaastra et al. 1999). Recently, XMM appears to have 
confirmed the presence of very extended, EUV/soft-X-ray 
excess emission in Abell 1795 (Arnaud 2000). In addition 
to the controversies about individual clusters, Arabadjis & 
Brcgman (1999) argue that all the detections may be due 
to uncertainties in Galactic absorption. 

Thus, it is not certain at this point whether excess EUV 
emission is a general characteristic of clusters of galaxies. 
In order to limit the contribution of cluster EUV emission 
to the extragalactic background, we will assume that it is, 
in fact, a general property of rich clusters. If this is the 
case, then our calculations may give a determination of 



the cluster contribution to the EUV background. If EUV 
excess emission is present in only a fraction of rich clus- 
ters, then our calculations will give a strong upper limit 
to the contribution to the extragalactic background. 

A promising theory for the EUV emission by clusters is 
that it is produced by inverse Compton (IC) scattering of 
the Cosmic Microwave Background (CMB) by relativistic 
electrons in the intracluster medium (Hwang 1997; Bowyer 
& Berghofer 1998; Enfilin & Biermann 1998; Sarazin & 
Lieu 1998). This requires the presence of a distinct pop- 
ulation of relatively low energy relativistic electrons with 
energies of ^100 McV (Bowyer & Berghofer 1998) beyond 
the population of higher energy electrons which produce 
radio emission in a small fraction of clusters (e.g., Feretti 
1999). Such a distinct population is expected theoretically, 
because the lower energy electrons have lifetime which are 
comparable to the Hubble time, while the radio emitting 
electrons have very short lifetimes (Sarazin 1999; Atoyan 
& Volk 2000; Takizawa & Naito 2000). 

An alternative theory is that the EUV emission is ther- 
mal in origin (Mittaz et al. 1998; Lieu, Bonamente, & 
Mittaz 1999a). There are a number of concerns with this 
thermal model, particularly with the energetics (Fabian 
1996; Sarazin & Lieu 1998). Recently, Lieu et al. (1999a) 
have argued that the observed EUV to soft X-ray spectra 
of clusters decline rapidly with increasing frequency, and 
that this favors a thermal model. However, such a rapid 
decline is, in fact, what was predicted by the IC model 
(Sarazin & Lieu 1998; Sarazin 1999; Atoyan & Volk 2000; 
Takizawa & Naito 2000; Figure 3 below). 

In any case, we will normalize our calculations of the 
EUV emission from clusters to the observed EUV flux from 
the Coma cluster, which is the best observed rich cluster. 
Thus, our results are somewhat independent of the emis- 
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sion mechanism for the EUV emission. In detail, we need 
a model to scale the EUV luminosity of Coma to other 
clusters and other redshifts, and we need a model for the 
spectrum of the emission. At present, there are no de- 
tailed thermal models for the EUV emission from clusters 
which provide such predictions. On the other hand, the 
IC model provides a very simple scaling of the luminosity 
and spectrum to other masses of clusters and to other red- 
shifts (§ 3 below). Thus, we will adopt the IC model for 
the emission to scale the observed flux of the Coma cluster 
to other clusters and redshifts. 

If this radiation really is a common feature of clusters at 
the indicated luminosity levels, it might make a contribu- 
tion to the intergalactic background radiation field at UV 
through soft X-ray energies. In this paper, we estimate 
the mean intensity of the diffuse EUV background from 
galaxy clusters, both at the current epoch and at higher 
redshifts. We will compare the radiation field from clusters 
with that produced by quasars and other active galactic 
nuclei (AGN) (e.g., Haardt & Madau 1996). In § 2, we 
use Press-Schechter (1974) theory to calculate the abun- 
dance of clusters of a given mass as a function of redshift. 
In § 3, the EUV emission of clusters of a given mass is 
calculated, based on a simple IC theory. We assume that 
both the thermal energy of the intracluster gas and the 
energy of relativistic electrons are produced ultimately by 
intracluster shocks, and that these shocks convert a fixed 
fraction of the shock energy into relativistic particles. The 
resulting diffuse radiation field at the present epoch and at 
higher redshifts are given in § 4. The resulting ionization 
rates of H I and He II are also derived there. The cluster 
diffuse EUV is compared to that produced by quasars and 
other AGN. The conclusions are summarized in 8 5. 



2. MASS DISTRIBUTION OF CLUSTERS 

2.1. Predicted Mass Distribution of Clusters 

In order to determine the total background radia- 
tion from galaxy clusters, it is first necessary to model 
their number density and its evolution with redshift z. 
While the present-day abundance of rich clusters is well- 
determined from observations (see § 2.2 below), the evo- 
lution of the cluster abundance is poorly determined at 
present, particularly at high redshifts. Thus, we will use 
the Press-Schechter (PS) formalism (Press & Schechter 
1974) to model the cluster density and its evolution, with 
parameters chosen to fit the observed present-day proper- 
ties of clusters and the results of more detailed numerical 
models of large-scale structure. Comparisons to observa- 
tions of clusters and to numerical simulations show that 
the PS provides an excellent representation of the statisti- 
cal properties of clusters, if the PS parameters are carefully 
selected (e.g., Lacey & Cole 1993; Bryan & Norman 1998; 
Viana & Liddle 1999). Let n(M,z)dM be the comoving 
number density of clusters with masses in the range M to 
M + dM in the Universe at a redshift of z. According to 
PS, this differential number density is given by 
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where ~p is the current mean density of the Universe, a{M) 
is the current rms density fluctuation within a sphere of 
mean mass M, and S c (t) is the critical linear overdensity 
for a region to collapse. In using this expression, we have 
made the usual assumptions that each cluster under con- 
sideration has recently collapsed (z b s « z co u) and that 
collapse is spherically symmetric. However, comparisons 
to numerical simulations show that this expression gives 
a good fit to the detailed results if the parameters are se- 
lected carefully (e.g., Lacey & Cole 1993; Viana & Liddle 
1999). We will assume a power-law spectrum of density 
perturbations, which is consistent with CDM models: 
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where as is the present day rms density fluctuation on 
a scale of 8 h- 1 Mpc, M 8 = (4tt/3)(8 h' 1 Mpc) 3 p is the 
mass contained in a sphere of radius 8 hr 1 Mpc, and the 
Hubble constant is H — 100 h km s _1 Mpc -1 . The expo- 
nent a is given by a = where the power spectrum of 
fluctuations varies with wavenumber k as fc n ; we assume 
n = —7/5, which leads to a = 4/15 throughout. Using the 
above relations, we can rewrite equation (1) as 
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The normalization of the power spectrum and overall 
abundance of clusters is set by as ■ Although there are pub- 
lished values for this constant (e.g., Viana & Liddle 1999; 
Bahcall & Fan 1998), we choose to normalize our calcu- 
lations to a fixed value for the observed local abundance 
of clusters. Specifically, for each cosmology, we choose 
the value of a s such that the present day integrated num- 
ber density of clusters with mass M > 8 x 10 14 h~ x Mq 
matches the observed value of 2 x 10~ 7 /i 3 Mpc~ 3 (Bah- 
call & Fan 1998). This method has the advantage that it 
forces the present day cluster abundance to be the same 
for each cosmology, as opposed to having a small spread in 
values, which occurs if an analytic approximation for the 
variation of as with cosmology is used. 

The cosmological model is characterized by several pa- 
rameters. The dependence on the Hubble constant H a = 
100 h km s _1 Mpc -1 just produces an overall scaling of 
the results, which we include through the factor h. We 
characterize the cosmological solution by the two stan- 
dard dimensionless parameters. First, fl = p/p c is the 
ratio of the current mass density to the critical mass den- 
sity p c = 3H$/(8wG). Second, fl A = A/(3H$), where A is 
the cosmological constant. We will consider three distinct 
cosmological models in our analysis: an Einstein-deSitter 
universe (fio = 1, = 0), an open universe (Clo = 0.3, 
Oa = 0), and a low-density flat universe (f2o = 0.3, 
Q\ = 0.7). The number density of clusters evolves dif- 
ferently in each of these three models. While the latter 
two models agree better with observations (e.g., Bahcall 
& Fan 1998), we include the closed model for complete- 
ness. 

The evolution of the density of clusters is encapsulated 
in the S c (z) term in equation (3). In general, 8 c (z) oc D(t), 
where D(t) is the growth factor of linear perturbations as 
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a function of cosmic time t (see Peebles [1980], § 11 for de- 
tails). Expressions for the 5 c (z) in different cosmological 
models are: 
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for the open, closed, and flat models respectively. For the 

— i ' — - 

open model (fl < 1, A = 0), t n = nH L Sl (1 - Q ) 2 
represents the epoch at which a nearly constant expan- 
sion takes over and no new clustering can occur, and the 
growth factor can be expressed as 

3 sinh r\ (sinh r) — rj) 



D(t) 
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where r\ is the standard parameter in the cosmic expansion 
equations (Peebles 1980, eqn. 13.10) 
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The solution for 5 C in the Einstein-deSitter model can 
be obtained from the open model solution by the limit 
tn/t — > oo (Lacey & Cole 1993). To evaluate 5 C in the 
flat model (Oo + = 1), we have used an approximation 
given by Kitayama & Suto (1996). Here 0/ is the value of 
the mass density ratio f2 at the redshift of formation, 
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In this model the growth factor can be written as 

D{x) = + 3/ 2 2 )1/2 [ X x 3 / 2 (x 3 + 2)- 3 / 2 dx (8) 

(Peebles 1980, eqn. 13.6) where x = (^f-) 1/3 and x = 
x /(l + z). 

Figure 1 shows how the cluster abundance evolves with 
redshift for the three cosmologies under consideration. For 
each model we plot the integrated number density of clus- 
ters at each redshift, which is simply the integral of equa- 
tion (3) over all masses M > 8 x 10 14 h^ 1 M©. The mass 
plotted here is the virial mass, which is the mass within the 
virial radius. The closed model shows the strongest evo- 
lution, followed by the flat model, with the open model 
showing the weakest evolution, as expected (e.g., Bahcall 
& Fan 1998). 

2.2. Comparison to Observed X-ray Luminosity Function 

As a check on these models, we have estimated the re- 
sulting X-ray luminosity function of clusters at the present 
epoch. To convert the predicted mass function of clus- 
ters to an X-ray luminosity function, we first determine 
the temperatures of the gas in clusters using the ob- 
served mass-temperature relationship (Evrard, Metzler, & 
Navarro 1996; Horner, Mushotzky, & Scharf 1999) 



M = c xl0 13 / l - 1 
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where T m = 1 keV, cq « 5, and c\ w 1.5. Then, the X-ray 
luminosities Lx were calculated from the X-ray tempera- 
tures using the the X-ray luminosity-temperature relation- 
ship (Arnaud & Evrard 1999) 



L x = b x 10 44 h- z [ ^ 



ergs 
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where Tx is the X-ray temperature, 7} = 6 keV is a char- 
acteristic value used for scaling, do « 2.88, and d\ w 2.88. 
Here, Lx is the bolomctric X-ray luminosity. Figure 2 
shows the resulting differential luminosity function at the 
current epoch as predicted by PS for our three cosmolo- 
gies. 

For comparison, we overplot the observed local bolomet- 
ric X-ray luminosity function from Ebcling ct al. (1997), 
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Here, L* x w 9.3 x 10 44 erg s _1 is a characteristic cluster 
X-ray luminosity, ao « 2.64, and a\ w 1.84. The model is 
in reasonable agreement with the observations. 



3. PREDICTED EMISSION FROM CLUSTERS OF GALAXIES 

As discussed in § 1, EUV and very soft X-ray emission 
has been detected from a number of clusters of galaxies. 
Although it is not yet certain how common this emission 
is, we will assume that it is a general feature of clusters. 
As we noted in § 1, this implies that our results are upper 
limits on the extragalactic EUV background due to clus- 
ters, if EUV emission is not a common phenomena. While 
the mechanism for this emission is still under debate, in- 
verse Compton (IC) scattering of the cosmic microwave 
background (CMB) by cosmic-ray electrons seems to be a 
very strong candidate (Hwang 1997; Bowyer & Bcrghofer 
1998; Enfilin & Biermann 1998; Sarazin & Lieu 1998). 
For the purposes of this analysis, we will assume that IC 
scattering is the dominant source of EUV radiation from 
clusters. However, we will normalize our models to the ob- 
served EUV flux of the Coma cluster; the IC model will be 
used to scale this flux to other cluster masses and rcdshifts. 

The IC luminosity from an individual cluster is given by 
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where Ucmb is the energy density in the CMB, Eqr is 
the total energy in cosmic-ray electrons in the cluster, and 
(7) is the average Lorentz factor of the cosmic-rays: 
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where ncR( r y)d , y is the number density of cosmic-ray elec- 
trons with Lorentz factors between 7 and 7 + dj. The 
EUVE observations of clusters suggest that (7) w 300 (for 
a Hubble constant of 50 km s _1 Mpc -1 ; Sarazin & Lieu 
1998). One expects that this value will correspond to elec- 
trons whose lifetimes are similar to the age of the cluster 
(Sarazin 1999), which implies that (7) oc h. Thus, we 
assume that 



(7) = 600/7 



(14) 



(Sarazin & Lieu 1998). The energy density in the CMB is 
well determined by observations, Ucmb — ^ cf tT cmb /c = 
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4.20 x 10~ 13 (1 + z) 4 ergs cm~ 3 , where T C mb is the tem- 
perature of the CMB, T C mb = 2.73(1 + z) K. 

The observed EUV luminosities of clusters imply that 
the total energies of cosmic ray electrons are Eqr ~ 10 62 
ergs (for a Hubble constant of 50 km s _1 Mpc ; Sarazin 
& Lieu 1998). These values are typically 1-10% of the 
thermal energy content of the intracluster gas, E gas . The 
observed values of E gas scale as E gas oc h~ 5 / 2 . The ther- 
mal energy content of the intracluster gas is the result 
of shock heating, where the shocks might include cluster 
merger shocks, overall accretion shocks, and shocks pro- 
duced by galaxy winds (Sarazin 2000). These shocks all 
have velocities of ^1000 km s _1 . Supernova remnants and 
other diffuse astrophysical shocks with similar velocities 
accelerate relativistic particles, and the radio properties of 
supernova remnants imply that ~3% of the shock energy 
goes into accelerating relativistic electrons (e.g., Bland- 
ford & Eichler 1987). Since both the relativistic electrons 
and the thermal energy in the intracluster gas arise from 
the same shocks, it seems reasonable to assume that the 
energy which is injected in relativistic electrons is propor- 
tional to the thermal energy content of the ICM. Thus, we 
will assume that 

E C R = fcREgas ■ (15) 
We adopt a fixed value for fcR which is determined by 
the observed EUV flux from the Coma cluster, which is 
the best observed rich cluster. To determine fcR, we ap- 
ply the Galactic absorption and redshift appropriate to 
Coma to our adopted spectral model (see Figure 3 be- 
low). Then, the resulting spectrum was convolved with 
the EUVE effective area to give the expected count rate. 
This was compared to the observed values of the EUVE 
count rate (Lieu et al. 1996a; Bowyer & Berghofer 1998; 
Bowycr ct al. 1999) which range over a factor of about 
two. The model spectrum was renormalized to give the 
observed count rate, and the total cosmic ray energy Eqr, 
was determined. The ROSAT X-ray observations of Coma 
were used to determine the total gas mass within the pro- 
jected radius of the EUV observations (e.g., Mohr, Math- 
iesen, & Evrard 1999). The corresponding thermal energy 
content of the gas E gas was derived. Finally, Eqr and 
Egas for Coma were compared to give fcR (eq. 15). The 
values were fcR, = (0.014 — 0.031)/i~ 1,/2 , corresponding to 
the range in observed count rates for Coma. We adopt the 
value fcR — 0.02 h^ 1 / 2 . This is consistent with the re- 
ported fluxes of other clusters (Sarazin & Lieu 1998). For 
any individual cluster, it is likely that Ecr and Jcr will 
depend on its dynamical history, but we assume that on 
average there is a simple relationship between the thermal 
energy in the gas and the energy in relativistic particles. 
The ICM thermal energy is given by 



E — 

^gas — 



3 kT 
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where T is the average temperature of the ICM, /i = 0.61 
is the mean mass per particle in units of the proton mass 
m p , and M gas is the mass of the ICM. There are a number 
of arguments which suggest that a portion of the thermal 
energy of the ICM may have come from non-gravitational 
effects which probably occurred before the clusters formed 
(e.g., Kaiser 1991; Cavaliere, Menci, & Tozzi 1997). Even 
if this heating involved shocks, it is unlikely that primary 



electrons accelerated in such shocks survive to the present 
epoch (Sarazin 1999). Thus, we will determine the gas 
temperature in our model clusters using a theoretical re- 
lationship based on purely gravitational effects, which is 
consistent with the PS formulation and with numerical hy- 
drodynamical simulations of cluster formation from large 
scale structure. We choose the mass-temperature relation- 
ship 

x [h 2 A c E(z) 2 ] 1/3 keV, (17) 

where fa ~ 0.8 is a normalization factor (Bryan & Nor- 
man 1998). The quantity E(z) is defined as E(z) 2 = 
Q {1 + zf + n R (l + z) 2 + fi A , where n R = l/(H R) 2 
where R is the current radius of curvature of space. A c 
represents the mean density of the cluster divided by the 
critical density at a given redshift. Fits for this parame- 
ter are given by Bryan & Norman (1998) for the relevant 
cosmological models: 



(18) 



18tt 2 + 82x - 39x 2 (Qr = 0) 
18tt 2 + 60x - 32x 2 (fi A = 0) 

where x = [Q (l + z f I E i z ) 2 ] - 1 - Note that Q R = is 
a consequence of the flat model (flat space has an infinite 
curvature). When Qq = 1, x = 0, giving A c = 187r 2 . 

Finally, we need to determine the total gas mass M gas 
of our clusters. X-ray observations indicate that the gas 
fraction in cluster is fairly constant for rich clusters, with a 
value of «22% (for a Hubble constant of 50 km s~ x Mpc -1 ; 
e.g., Arnaud & Evrard 1999). Thus, we assume that 

M gas =f gas M, (19) 

with f gas = 0.07 h- 3 / 2 . 

Using equations (12)-(19), we calculate the average IC 
luminosity Lie of a cluster of a given total mass M. For 
individual clusters, one expects the spectrum of IC EUV 
emission will depend on the history of particle acceleration 
in the cluster (Sarazin 1999). For this study, we adopt a 
"typical" IC EUV spectrum, since we intend to average 
over an ensemble of clusters to determine the cluster con- 
tribution to the EUV background. Specifically, we assume 
the spectral shape given for Model 11 in Sarazin (1999). 
This adopted model spectrum L u (model) is shown in Fig- 
ure 3, where L v is the luminosity per unit frequency v. 
This particular model spectrum has a total IC luminosity 
of Ljc'(model) = 2.86 x 10 44 ergs s _1 . The normalization 
of the spectrum of each model cluster is scaled to give the 
correct value of Lie for that cluster. At redshifts z > 0, 
the spectrum is blue shifted due to the increase in the tem- 
perature of the CMB. Thus, the spectrum of emission of 
a cluster with total mass M at redshift z is given by 



L IC (M,z,u) 



L v/[x+z) (model) 

Lic{M,z)— — — — . (20) 

(1 + z)Lic(model) 



3.1. The Diffuse Radiation Field 

The total emissivity of cluster IC emission at any given 
redshift is found by integrating the emission due to clus- 
ters of a given total mass M (eq. 20) over the PS mass 
function (cq. 3): 

e{v,z) = j L IC (M,z,v)n{M,z)dM . (21) 
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The total emissivity over all frequencies is just 

e t ot(z) = J e{v,z)dv = j L IC (M, z) n(M, z) dM . (22) 

The equation of radiative transfer for the emission gives 
a mean intensity for the radiation field at an observed red- 
shift z a and observed frequency v Q of (Haardt & Madau 
1996) 

dl (1 + z f 



1 f°° 
J{v ,z ) = — J dz- 



dz (1 + z) 3 

xe(i/, z)e- T «"(""' z ">*>. (23) 

Here, dl / dz is the line element in a Friedmann cosmology, 
given by 
dl_ 

dz 



H 



x [n (i + zf 

The quantity T e ff(v ,z n 



f n fl (i + zf + n A ] 1/2 . (24) 

z) gives the effective optical 
depth of the intergalactic medium (IGM) at the observed 
frequency between z and z. We can simplify equation (23) 
by noting that the energy density in the CMB, which is 
present in our expression for e(z) through the Ljc term, 



goes as (1 + z) /(l 



When we consider the emis- 



sivity per unit frequency e(u, z), we get another factor of 
(1 + z )/(l + z) due to the redshift applied at each fre- 
quency. Multiplying these factors together gives us a net 
factor of (1 + z) 3 /(l + z Q ) 3 , which directly cancels the red- 
shift term in equation (23). We can therefore rewrite this 
equation as 



J{v , z ) 



1 

4"7T 



dz ^- e(u , z) e 
dz 



f{v ,Z ,z) 



(25) 



where e{v a , z) is defined as the emissivity for a cluster pop- 
ulation at a redshift z calculated as though it were actually 
located at redshift z Q . 

The opacity of the IGM is mainly due to the effective 
column density of the less ionized clouds (such as Lya 
clouds). In determining the optical depth of the IGM, we 
closely follow the treatment in Haardt & Madau (1996). 
They write the optical depth of the IGM as 

d 2 N 



T e ff{Vo, 



z , z)= j dz' / 

Jz„ JO 



dV 



HI 



l-e 



-t(N hi ,z' 



dN HI dz' 



(26) 



where Nhi is the column density of neutral hydrogen in 
a IGM cloud, and {d 2 N / dN HI dz') gives the number of 
clouds per unit H I column and per unit redshift z' . The 
quantity t(Nhi , z' , v) is the optical depth of a cloud with 
an H I column of N^i at redshift z' for the frequency 
v = u {l + z')j{\ + z Q ); this depends on the helium abun- 
dance and the ionization structure of the clouds; we adopt 
the abundances and ionization structures calculated in 
Haardt & Madau (1996). For the ionization, we use the ap- 
proximation given in their equation (12), with rjthin taken 
from their Figure 7. 

4. RESULTS 

4.1. Models without IGM Absorption 

We first calculate the diffuse radiation field neglecting 
IGM absorption (r e // = in cq. [23]). Here, we consider 



only the values at the present epoch, z = 0. We find the 
total mean intensity to be 1.0 x 10~ 10 ergs cm~ 2 s _1 sr _1 
for the open model, 9.0 x 10 -11 ergs cm~ 2 s _1 sr _1 for 
the flat model, and 4.7 x 10~ n ergs cm~ 2 s _1 sr _1 for the 
closed model. 

We have calculated the spectrum of the cluster EUV 
background at the present epoch, z = 0, neglecting absorp- 
tion for our three cosmological models. Results are plotted 
in Figure 4. The corresponding ionization rates have been 
determined by integrating each of these radiation fields 
over the ionization cross-sections of H I and He II. The 
results are ( H = 1-8 x 10~ 17 , 1.6 x lO" 17 , 8.6 x 10~ 18 and 
ChcII = 2.05 x 10~ 18 , 1.8 x 10~ 18 , 9.7 x 10~ 19 s" 1 for the 
open, flat, and closed cosmological models, respectively. 



4.2. Models with IGM Absorption 

Next, we calculate models including the opacity of the 
IGM, following the treatment of Haardt & Madau (1996). 

4.2.1. Present Epoch Diffuse Radiation 

We have first calculated the values for the diffuse clus- 
ter IC radiation field at the present epoch, z = 0. We 
find the total mean intensity to be 9.0 x 10 -11 ergs cm~ 2 
s _1 sr _1 for the open model, 8.2 x 10~ n ergs cm" 2 s" 1 
sr _1 for the flat model, and 4.5 x 10~ n ergs cm~ 2 s _1 
sr _1 for the closed model. We note that the numbers 
are only slightly smaller than the values ignoring absorp- 
tion (§ 4.1 above). Unlike the ionizing background from 
quasars, most of the diffuse emission from clusters is pro- 
duced at relatively small redshifts (§ 4.2.2 below), where 
the IGM opacity is low. 

The predicted spectrum of the EUV background from 
clusters at the present epoch z Q = is shown in Fig- 
ure 5 for each of our three cosmological models. The 
"bumps" in these curves are due to absorption by H I 
and He II. Clearly, absorption has only a small effect on 
the spectrum (compare with Figure 4). The present day 
ionization rates for cluster EUV models with IGM absorp- 
tion arc ( H = 1.6 x 10~ 17 , 1.5 x 10" 17 , 8.2 x 10~ 18 and 
( HeII = 1.6 x 10~ 18 , 1.4 x 10~ 18 , 8.5 x 10- 19 s" 1 for the 
open, flat, and closed cosmological models, respectively. 
Ionization rates from QSOs are given below for compari- 
son (§ 4.2.2). 

Figure 6 compares the predicted contributions of clus- 
ters and quasars to the diffuse EUV background at z D = 
for the open cosmological model. The quasar results are 
taken from Haardt & Madau (1996), Figure 5a. A compar- 
ison of the two curves shows that the contribution to the 
EUV background from clusters is < \% of that of quasars. 
Thus, while clusters do make a contribution to the EUV 
background at the present epoch, the predominant source 
is due to AGN. 

4.2.2. Redshift Evolution 

The spectrum of the diffuse radiation field due to clus- 
ter IC EUV is plotted as observed at redshifts of z a = 0, 
0.5, and 1.0 for the open cosmological model in Figure 7. 
Models were calculated for higher redshifts {z > 1.5), but 
the values much smaller than those shown. This figure 
clearly shows that the intensity of the EUV background 
from clusters decreases strongly with increasing redshift, 
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which means that the largest contribution to the EUV 
background from clusters occurs at the current epoch. On 
the other hand, the predicted intensity of the diffuse UV 
radiation field due to quasars increases with increasing red- 
shifts out to z a ~ 3. (Haardt & Madau 1996, Figure 
5). Thus, the fractional contribution to the EUV back- 
ground from clusters as compared to quasars also decreases 
strongly with redshift. This is because the abundance of 
rich clusters increases with time due to gravitational merg- 
ers of smaller systems, while quasars were most abundant 
in the early universe (at z > 2). 

The ionization rate due to cluster EUV also decreases 
rapidly with increasing redshift. We find rates of = 
1.6 x 1(T 17 , 2.0 x 1CT 18 , 1.2 x 1(T 19 , 3.9 x 1CT 21 s" 1 
for H I and C, HeII = 1.6 x 1(T 18 , 2.1 x 1(T 19 , 1.3 x 1CT 20 , 
4.7x 1(T 22 s" 1 for He II at z a = 0, 0.5, 1.0, and 1.5, respec- 
tively, for the open cosmological model. To compare these 
results with those from QSOs, we use the results given by 
Haardt & Madau (1996, Fig. 6) for an open cosmologi- 
cal model with qo = 0.1, which give (h = 4.1 x 10~ 14 , 
1.6 x 10~ 13 , 4.5 x 10~ 13 , 9.3 x 10~ 13 s" 1 for H I and 
(Hen = 4.1 x 10~ 16 , 1.6 x 10" 15 , 4.4 x 10" 15 , 9.0 x 10" 15 
s _1 for He II at z Q = 0, 0.5, 1.0, and 1.5, respectively. In 
general, the ionization rates due to cluster EUV are <1% 
of those due to QSOs. The largest contribution is for the 
He II ionization rate at z = 0, where the cluster EUV gives 
^0.4% of the quasar rate. 

5. CONCLUSION 

Observations with EUVE indicate that some clusters of 
galaxies are luminous sources of extreme ultraviolet (EUV) 
radiation. In order to limit the contribution of EUV emis- 



sion from clusters to the extragalactic background, we have 
assumed that this emission is a common feature of clus- 
ters. A promising theory for this emission is that it is 
due to IC scattering of the Cosmic Microwave Background 
photons by relativistic electrons. We have given a sim- 
ple model for the average EUV luminosity expected from 
clusters based on this theory. We summed this emission 
over clusters of varying mass at different redshifts, using 
the Press-Schechter formalism to determine the mass func- 
tion of clusters as a function of redshift. We determined 
the amount of background radiation produced by clusters 
through EUV IC emission. The total mean intensity, spec- 
trum, and the ionization rates for H I and He II were 
determined at present and at a variety of redshifts. Be- 
cause clusters form by the merger of smaller subclusters 
and the abundance of massive clusters increases with time, 
the amount of EUV background radiation should be larger 
at present than in the past. We compared our results to 
the ionizing background expected from quasars. We find 
that while clusters do contribute a significant EUV back- 
ground, it is less than a percent of that expected from 
quasars. Of course, this last conclusion is strengthened if 
EUV emission is not a general feature of clusters. 
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Fig. 1. — The evolution of the abundance of clusters with redshift z for three distinct cosmologies. Only clusters with mass M > 
8 X 10 14 h~ 1 Mq within their virial radius R v i r are considered. The solid curve represents an open Universe (Qq = 0.3 , Q\ = 0), the dashed 
curve a flat Universe (Qq = 0.3 , S1a = 0.7), and the dotted curve a closed Universe (Qg = 1.0 , SIa = 0). 
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L x (h" 2 1 44 ergs/sec) 

Fig. 2. — The present day cluster bolometric X-ray luminosity function. The solid curve is the observed luminosity function (Ebcling ct 
al. 1997). The dashed curve shows the predicted luminosity function for a closed Universe (f2o = 1, Qa = 0), while the dotted curve shows 
the predicted luminosity function for an open Universe (f!o = 0.3, SI a = 0). The flat Universe is omitted here for clarity since the predicted 
curve essentially matches the predicted curve for the open Universe. 
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Fig. 3. — The adopted IC EUV spectrum for our clusters, which is the spectrum of Model 11 in Sarazin (1999). Here, L„ is the luminosity 
per unit frequency v. This model has a total IC luminosity of Ljq = 2.86 X 10 44 ergs s — 1 . The normalization of the spectrum of each model 
cluster is scaled to give the correct value of Ljc for that cluster. 
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Fig. 4. — The mean specific intensity of the EUV background expected from clusters at the current epoch for three distinct cosmologies. 
The solid curve represents an open Universe (do = 0.3, SIa = 0), the dashed curve a flat Universe (Ho = 0.3, SIa = 0.7), and the dotted 
curve a closed Universe (Qq = 10, Qa = 0). Absorption effects are ignored. 
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Fig. 5. — The mean specific intensity of the EUV background predicted by models including IGM absorption by H I and He II. The notation 
is the same as in Figure 4. 




Fig. 6. — The predicted mean specific intensity of the EUV background at z = for an open Universe. The curve labeled "Clusters" shows 
the contribution from clusters while the curve labeled "AGN" shows the contribution from quasars. 
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Fig. 7. — The spectrum of the diffuse EUV background from clusters as observed at rcdshifts z = 0, z = 0.5, and z = 1 for the open 
cosmological model including IGM absorption. 



